;$Id: lognormal_r.pro 4 2009-07-01 08:10:55Z dingyy.pku $
FUNCTION LogNormal_r, N, d_bar, d, sigma
  pi=3.14159265D;
  logSigma=Alog(sigma)
  twoPiSqrt=Sqrt(2.0d*pi);
  result=N/(twoPiSqrt*logSigma*d)*Exp(-(Alog(d)-Alog(d_bar))^2/(2.0d*logSigma^2));
  Return, result;
END
